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Введение. Рассмотрен тепловой расчет объемной конструкции с использованием метода конечных элементов. 
По планам Министерства энергетики РФ в стране в ближайшие годы будет создана мощная индустрия ветровой 
энергетики. В связи с этим расчеты при производстве строительных конструкций ветроэнергетических устано- 
вок в настоящее время приобретают большое значение. Производство таких конструкций из стеклопластика 
представляет собой сложный термохимический процесс, включающий полимеризацию связующего при строго 
заданных тепловых режимах. Целью работы является разработка методики трехмерного конечно-элементного 
расчета нестационарного режима нагрева составной конструкции сложной формы. 

Материалы и методы. Определение полей температур конструкции сложной формы, выполненной из неодно- 
родных материалов, приводит к необходимости применения численных методов и, в первую очередь, метода 
конечных элементов. Конечно-элементное моделирование поведения композиционных материалов при формо- 
вании до настоящего времени остается незавершенным. Для частичного ее решения выполнена адаптация из- 
вестного уравнения теплопроводности для конкретной задачи на основе первого закона термодинамики. Пред- 
ложены новые конечно-элементные модели, описывающие тепловые поля в конструкции при ее изготовлении. 
Определена точность моделирования тепловых процессов. Проведено численное моделирование нагрева. 
Результаты исследования. Решение задачи выполнено в многофункциональном программном комплексе 
АМЬУБ с реализацией методики расчета на языке параметрического программирования АРПОГ.. Рассчитаны по- 
ля температуры лопастных элементов ветроэнергетических установок на этапе их изготовления, что позволило 
выявить характерные особенности технологического процесса производства данных конструкций и получить 
рекомендации по уточнению процесса их склейки. 

Обсуждение и заключения. Полученные результаты могут быть использованы в тепловых расчетах элементов 
строительных слоистых конструкций сложной формы из композиционных материалов в ветроэнергетике, ма- 
шиностроении, авиастроении, судостроении, приборостроении ит. д. 


Ключевые слова: конечно-элементный расчет, температурное поле, неметаллическая конструкция, технологи- 
ческий процесс, моделирование. 
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Введение. Тепловые расчеты при производстве строительных конструкций ветроэнергетических уста- 
новок в настоящее время приобретают большое значение. По планам Министерства энергетики РФ' в стране в 
ближайшие годы будет создана индустрия ветровой энергетики. В 2016—2017 гг. в нее пришли крупные россий- 
ские и иностранные инвесторы, которые взяли обязательства по развитию технологической и производственной 
базы. Устраняются основные недостатки лопастных ветроэнергетических установок путём совершенствования 
их конструкции и технологии изготовления. В первую очередь это относится к лопастям, которые являются 
основными элементами ветротурбины, концентрирующими главную интеллектуальную составляющую уста- 
новки. Производство стеклопластиковой лопасти представляет собой сложный термохимический процесс, 
включающий полимеризацию связующего при строго заданном тепловом режиме. В связи с этим точный рас- 
чет температурных полей в теле лопасти при ее формовании имеет большое практическое значение. 

Определение полей температур конструкции сложной формы, выполненной из неоднородных материа- 
лов, приводит к необходимости применения в расчетах численных методов и, в первую очередь, метода конеч- 
ных элементов (МКЭ). Проблема конечно-элементного моделирования композиционных материалов при фор- 
мовании остается до настоящего времени нерешенной. Это определяется сложностью поведения композита как 
неоднородного материала, обладающего ярко выраженными термозависимыми физическими, в т. ч. реологиче- 
скими свойствами, остаточными внутренними явлениями, старением и др. 

Современные модели учитывают анизотропность материала, его пластическое поведение при сложном 
напряженно-деформированном состоянии и др. Так, в работах [1, 2] исследованы температурные свойства ком- 
позиционных материалов. Для обоснования и подтверждения влияния распределения и ориентации упрочняю- 
щих частиц на коэффициент температурного расширения материала авторами выполнено конечно-элементное 
моделирование. При геометрическом построении за основу была принята реальная микроструктура. Расчетами 
установлено, что экспериментально обнаруженная анизотропия температурных свойств объясняется внутрен- 
ними напряжениями в композите, зависящими от ориентации наполнителя. 

Многие исследования посвящены расчетам двумерных моделей армированных композиционных мате- 
риалов. Для анализа и оптимизации свойств композита предложена расчётная конечно-элементная схема, учи- 
тывающая его температурные свойства [3]. С помощью МКЭ рассчитаны поля остаточных напряжений, обра- 
зующихся в матрице в результате термической обработки для модели композита, включающего армирующие 
волокна [4]. С использованием конечно-элементной модели рассчитаны эффективные упругие характеристики 
многослойного композиционного материала, причем каждый слой описывался собственными термомеханиче- 


'Возобновляемая энергетика в России состоялась как отрасль / Министерство энергетики российской федерации // шшепегоо.соу.та : [сайт]. 
ОВГ: 66рз://ллщшепегоо.соу.га/по4е/10307/ (дата обращения: 04.02.2021). 
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скими свойствами [5]. Моделирование поведения анизотропных материалов и другие вопросы термонагруже- 
ния проанализированы также в работах [6-13]. 

Большинство исследований, базирующихся на конечно-элементном анализе, выполнены с использова- 
нием универсальной программной системы АМЗУЪ, позволяющей решать линейные и нелинейные, стационар- 
ные и нестационарные пространственные задачи теплопередачи и теплообмена. Например, с помощью АМЗУЗ 
решены задачи: 

® нестационарного нагрева конструкции при наличии зависимости от времени коэффициента конвекции”; 

® теплового расширения элементов при интенсивном нагреве”: 

® тепловыделения при пластической деформации"; 

® передачи теплоты в многослойных конструкциях через границу раздела сред” и др. 

Исследованию процессов обработки полимерных композиционных материалов посвящены работы 
С. Н. Шевцова, например [14]. Несмотря на большое количество работ, посвященных численным тепловым 
расчетам, актуальность оценки теплового режима формования конструкций лопастных установок сохраняется. 
Повышение требований к качеству производства элементов ветротурбин приводит к необходимости дальней- 
шего уточнения математических конечно-элементных моделей теплопередачи на основе более полного учета 
внешних и внутренних факторов. 

Целью работы является исследование нестационарных полей температур составной конструкции слож- 
ной объемной формы лопастного элемента ветроэнергетической установки в процессе ее изготовления при учете 
физических особенностей и геометрических нелинейностей. Одна из конструкций состоит из лонжерона и хвосто- 
вого отсека, части которых выполнены из стеклопластика, резины, вспененного пластика и др. Хвостовой отсек 
изготавливается по клеевой технологии в специальном приспособлении, которое состоит из верхнего и нижнего 
полукорпусов с мощными ребрами жесткости. Между ними закладывается лопастной элемент. Для его качествен- 
ной склейки необходимо строго выдерживать скорость нагрева и температуру выдержки в клеевых швах. Задание 
теплового режима осуществляется размещением приспособления с отсеком в нагревательной печи. 

Постановка задач работы состоит в следующем. Требуется создать конечно-элементную модель нагре- 
ва лопастного элемента при его формовании, рассчитать объемные поля температуры в любой момент процес- 
са, определить точность моделирования тепловых процессов в теле конструкции, а также сформировать график 
набора температуры в печи, обеспечивающий требуемый режим нагрева клеевых зон отсека. 


Материалы и методы. Основой тепловых расчетов является утверждение первого закона термодина- 
мики о том, что количество теплоты, полученное изолированной системой, затрачивается на изменение ее внут- 
ренней энергии. Применительно к элементарному объему это утверждение можно выразить математически: 


Г ТГ 
рс(ОТ / б-+ {У} ИТ) {Г Ч =а,, (1) 
где р — плотность материала; с — удельная теплоемкость; Т — температура; Е — время; {а} — вектор теплово- 
го потока. Чу —_ тепловыделение единичного объема. 


Векторный оператор И, и вектор скорости {7 для теплопередачи, реализуемой массой: 


в 
а 
Г 


где х, у, 2 — координаты. 
По закону Фурье тепловой поток связан с градиентом температуры: 


ч=-[РАЕТ, (2) 
где [2] — матрица теплопроводности: 
° Треппа! Тиле-Тгапуеп Тоадшо ап боюбоп ш Апзуз / ЗипаТесНОгоар // зшииесКеточр.сот : [сайт]. — ОВГ: 
Брз://\у\илу. зпилиесВогопр.сот/ар$-ап-васК$/еа-а 1с]е$/97-Реа-Нрз-блск$-Фегта]-гапз1епе. (дата обращения: 04.02.2021). 
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Ехцете ТВепта| Ехрапзюп Модейи> ш Апзуз Месвашса[ У’ойфепсь / ЗипаТесвОтоир // зпищесвотоир.сот : [сайт]. — ВП: 
В рз://Л\у\миму. зшишесйетоир.сот/йрз-ап4-влск$/еа-агсе$/139-ехгете-@егтпа!-ехрапзюп-то4еНпе-т-апзуз-теспашса|. (дата обращения: 04.02.2021). 
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сепегайоп-р!а5ис-4еогтаНоп. (дата обращения: 04.02.2021). 
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[р]=ол, о 
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Объединяя (1) и (2), получим: 
(ОТ / +4 (ПТ) = М (ПРИЫТ)+а,. (3) 


Раскрывая (3), получим: 


ОТ ОТ ОТ ОТ д ОТ) д ОТ] д ОТ 
ре(— Г. — +И, — + —)=а +—| А — +— А, — |+—| А, — 
0: бх ‘^ ду 02 Ох Хх) д\ ’д] 42 02 
Основное уравнение теплопроводности дополняется граничными условиями: 
|. Задание температуры на поверхности: Т = 7*. 


2. Задание теплового потока на поверхности: 


Т * 
{9} п}=-9, 
где {и} — нормальный к поверхности единичный вектор; а* — удельный тепловой поток. 
3. Теплопередача соответствует закону Ньютона (поверхностная конвекция): 


{@} {п} =а, (Т,-Т,), 


где а/— коэффициент теплопередачи, Т, — температура на поверхности, 7» — температура граничащей среды. 
Для случая постоянных теплофизических коэффициентов исследуемые тепловые процессы описывают- 
ся трехмерным уравнением нестационарной теплопроводности (уравнением Фурье): 
ОТ СТ СТ ОТ 
01 0х д’ 02 











где Л, — коэффициент теплопроводности. 

В связи с симметричностью задачи расчет проводится лишь для половины конструкции. Граничные 
условия для плоскости сечения определяются из условия: 4 = 0. 

На стыковочных поверхностях предполагается идеальный тепловой контакт, отсутствие термического 
сопротивления: 


ОТ 


2 


РТ 








+0 > +0 ° 


ОГ 
^. о А 
Ох 

Конвективные тепловые потоки подводятся к наружным поверхностям приспособления. Внутреннее 
тепловыделение при нагреве не учитывается. Начальные условия предполагают фиксацию постоянной темпе- 
ратуры по всей исследуемой области: Т = Ту. 

Задача решена методом конечных элементов в многофункциональном программном комплексе АМЗУЗ 
на языке параметрического программирования АРОГ. Конструкция нагревается с наружной поверхности 
конвекцией. Результаты разбивки приспособления и отсека на блоки, а также на объемные тетрагональные 
конечные элементы представлены на рис. 1. 





а) 6) 6) 


Рис. 1. Склеечное приспособление с хвостовым отсеком: а — общий вид; 6 — поверхности модели с указанием их номеров; 


в — разбивка анализируемой области на конечные элементы 


Теплофизические свойства материалов, используемых в расчетах, приведены в таблице 1. 
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Таблица 1 


Теплофизические свойства материалов 


ВТК) хД т 
[236 гб 





798 57 


Созданная для решения поставленной задачи программа определяет значения температуры в каждом 
расчетном узле на протяжении всего технологического процесса склейки. 

Алгоритм расчета состоит из следующих шагов: 

® задание исходных данных (геометрии, теплофизических свойств, начальной температуры ит. п.); 

® разбивка области на конечные элементы; 

® приложение тепловых источников; 

® определение температуры области в конце шагов нагрева; 

® вывод результатов расчета (тепловых потоков, температурных полей) на печать. 

Для проверки правильности работы программа проходила тестирование. На рис.2 показаны 
экспериментальные и расчетные зависимости температуры от времени в клеевой зоне (нижние кривые) и в 
самой печи (верхние кривые). Абсолютная погрешность расчета температуры в точке склейки не превышает 
5 °С, относительная погрешность для реально значимых температур — не более 6 %. 
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Рис. 2. Зависимости от времени расчетных и измеренных температур в контрольной и задающей точках 


Результаты исследования. Технологический процесс включает в себя монотонный, непрерывный 
нагрев склеечного приспособления и лопастного элемента до определенной температуры и ее поддержание за- 
данное время. Результаты расчета приведены на рис. 3-7. Так, на рис. 3 представлено распределение по объему 
в определенный момент времени тепловых потоков в векторной форме и суммарное количество теплоты в теле 
приспособления с отсеком. 
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264 Рис. 3. Тепловые потоки (а) и суммарный нагрев (6) приспособления с отсеком 
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Так как скорость нагрева приспособления относительно невелика, тепловые потоки перераспределяют- 
ся по его объему равномерно и оказываются приблизительно одинаковыми по всем направлениям. Небольшие 
исключения касаются лишь ребер жесткости. Конвективные тепловые потоки подводятся с боков, имеющих 
большую площадь, а отводятся теплопроводностью в тело полукорпусов через узкое поперечное сечение. 

На рис. 4 представлены температурные градиенты в теле приспособления с отсеком в один из моментов 
нагрева. Видно, что тепло распространяется от наружных поверхностей приспособления к внутренним обла- 
стям. Тепловые потоки наиболее активно подводятся в местах близких к боковому вкладышу и перемещаются 
одновременно как от бокового вкладыша, так и от верхнего и нижнего полукорпусов приспособления. На этом 
же рисунке представлен объемный вид температурного поля в аналогичный момент времени. Поскольку ло- 
пастной отсек неметаллический, его теплопроводность низка, и распределение температуры оказывается не- 
равномерным. 
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а) 0) 
Рис. 4. Векторные поля градиента температуры (а) и температуры (6) в теле приспособления с лопастным отсеком 
Зависимости температуры от времени в отдельных точках клеевой зоны представлены на рис. 5—7. Из- 
менением значений задающих температур окружающей среды и количества их переключений можно добиться 


максимального соответствия графика набора и выдержки температуры в клеевых зонах требованиям техноло- 
гического процесса. 
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Рис. 5. Зависимости температуры от времени в характерных точках клеевых зон и в блоке (зеленая линия) для технологиче- = 
ского режима, включающего два переключения при задающих температурах 160 и 140 °С > 
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Рис. 6. Зависимости температуры от времени в характерных точках клеевых зон и в блоке (зеленая линия) для технологиче- 
ского режима, включающего четыре переключения при точности задания температуры до 10 °С 
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Рис. 7. Зависимости температуры от времени в характерных точках клеевых зон и в блоке (зеленая линия) для технологиче- 
ского режима, включающего три переключения 


Обсуждение и заключения. Выполненный тепловой расчет элементов строительных конструкций вет- 
роэнергетических установок на этапе изготовления позволил выявить характерные особенности нагрева и по- 
лучить рекомендации по уточнению параметров технологического процесса их склейки. 
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